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^ ■ We present a first study of gluonic twist-four corrections to the deep inelastic struc- 

j^ . ture function F2 at small x and small Q^. The calculations are based upon the double 

logarithmic approximation of the coupled twist-four evolution equations of the gluonic 
twist-four operators that are expected to be dominant at small x. We first review the 
analytical results which are presently available and define the framework of our calcu- 
lation. In the second part we discuss the connection with DIS diffractive dissociation 
which can be used to estimate the size of some of the twist-four corrections. In the 
final part we show, for three different choices of the input distributions, the relative 
magnitude of the leading- twist and the twist-four contributions. 
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1 Introduction 

The observation that the rise of the deep inelastic structure function F2 at smah x starts already 
at surprisingly low values of Q^, has stimulated attempts to use the DGLAP evolution equations 
in the low Q^ kinematical regime [^0, where previously - i.e. before the advent of HERA - 
perturbative QCD had not been expected to be applicable. As a general result of these studies, 
with a suitable choice of initial conditions (in particular: not rising at small x) it is possible to 
describe, within the framework of the leading- twist, next-to-leading-order QCD evolution equations, 
the Q^-evolution of the structure functions down to rather low Q^ values. However, recently Q 
evidence has been given that a rather sharp transition away from this DGLAP evolution occurs 
between 3 GeV^ < Q^ < SGeV^, 10"^ < x < 10"^ 

The apparent success of DGLAP in the low Q^ region, however, does not provide much help 
for answering the question of how the transition from the perturbative parton picture into the 
non-perturbative hadronic region at Q^ = works. Most naively, one expects the breakdown of 
perturbative QCD to be accompanied by the growth of non-leading perturbative terms. For deep 
inelastic scattering, the expansion parameters are powers of the strong coupling constant as and 
powers of m? /Q"^ (with m denoting some hadronic scale) which define the twist expansion. Whereas 
the former expansion has already received some attention (e.g. the question of resummation in the 
anomalous dimension for small x [^), the role of non-leading twist at small x so far has not yet 
been addressed. In particular, if at low Q^ the transition between the parton picture and Regge 
physics results from an interplay of leading- twist and higher- twist corrections, one should expect 
to see the presence of a negative twist-four term already slightly above the transition. The results 
of Q indicate that twist-four might be present even at not so small values of Q^. If this is the 
case, existing leading-twist parametrizations in the low Q^-region need to be re-examined. If, on the 
other hand, higher- twist really remains small down to unexpected low values of Q^, one should look 
for an explanation of this phenomenon. In any case, a closer investigation of the role of higher-twist 
in the small x, low Q^-region is quite important. 

Twist-four operators are known Q to have their own evolution equations. For the small x domain, 
gluon operators are expected to be the dominant ones. Whereas fermionic twist-four operators 
have been investigated in some detail, very little is known about the gluonic ones. Following 
the arguments of [^ we expect that twist-four gluonic contributions are obtained by investigating 
the Q^-dependence of QCD diagrams with two, three, or four gluons in the t-channel. Explicit 
calculations of these diagrams have been done only for the small x behavior of gluonic scattering 
amplitudes: the BFKL |8[ (ladder) amplitude, from which one obtains the high energy behavior 
of the elastic scattering of a virtual photon on a gluon, and diffractive deep inelastic scattering 
in the triple Regge region |^], which allows to extract contributions with three and four gluons 
in the t-channel. Expanding these amplitudes in inverse powers of Q^ we expect to reproduce 
the Q^ evolution equations for twist-two and twist-four in the double leading log approximation 
(DLA): leading logarithmic in both InQ^ and Inl/j;. A complete leading order calculation of the 
Q^-evolution equations of the twist-four gluon operators is missing, and it is very important that 
this task will be addressed. In this paper we will adopt the viewpoint that the combined limit: 
a; ^ 0, Q^ — > cxD of diagrams with two, three, or four gluons in the t-channel provides the DLA for 
the twist-four gluon operators. 

A numerical analysis of twist-four corrections to, say, F2 follows the same sequence of steps as the 
usual twist-two strategy: for an input scale Qq one has to choose an initial distribution (with a 
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certain number of free parameters), and by application of evolution equations one gets values of 
the structure functions for larger Q^. Since the number of twist-four operators that mix under 
renormalization will, in general, be larger than in the twist-two case, the number of free parameters 
in such an analysis will be large, and a reliable combined analysis of twist-two and twist-four terms 
looks extremely difficult. It might, therefore, be helpful to include information on specific final 
states, which are known to belong to higher-twist. At present the most promising candidate are 
diffractive final states, in particular the diffractive production of vector mesons from longitudinal 
photons: this process has been shown to be calculable within perturbative QCD |10-^]. It belongs 



to higher- twist, and the measured cross section is substantial (about 20 % of the total diffractive 
cross section). Another potential source of information are diffractive final states with hard jets. 

In this paper we perform a first numerical estimate of gluonic twist-four corrections at small x 
and low Q^. For this we make use of the currently available elements of the twist-four evolution 
equations. In our scheme we, first, review the existing small x calculations of gluonic scattering 
amplitudes, and extract from them the relevant pieces of the Q^ evolution. The currently existing 
calculations do not allow to go beyond DLA. Generally we are dealing with gluonic amplitudes 
which couple to the photon through a fermion-loop. As a first general result we find that, for 
a systematic study of twist-four, we have to consider both the transverse and the longitudinal 
polarization of the photon. The reason for this is the fact that in the fermion-loop the twist-four 
contribution of the transverse photon looses one power of In Q^ , whereas the longitudinal photon 
does not. This is quite opposite to the leading- twist case where the transverse photon results in 
a logarithm, and the longitudinal one does not. Unfortunately, in DLA this result forces us to 
compare contributions of different order in as'- in order to make a legitimate comparison between 
the twist-four transverse and longitudinal contributions, we should know the longitudinal one up to 
the (in the sense of DLA) next-to-leading terms. We, therefore, restrict ourselves to ratios of twist- 
four corrections to DLA leading- twist results. A second general result of our investigation is the 
sign structure of the different higher- twist terms. Both for the longitudinal and for the transverse 
twist-four contributions we shall argue that there are four different pieces that have to be taken 
into account, and they come with alternating signs. The total twist-four contribution, therefore, 
depends crucially on the magnitude of each individual term, i.e. the size of the initial condition. 

The second purpose of this paper is the search for a practical procedure for estimating the higher- 
twist contributions. We propose to use experimental data of DIS diffractive final states. It is, 
therefore, necessary to first investigate the role of higher-twist in diffractive dissociation and to 
analyse the connection with the inclusive structure functions. In the second part of this paper 
we, therefore, review diffractive dissociation, with particular emphasis on higher-twist. For the 
determination of free parameters in our numerical analysis we will make use only of longitudinal 
quark-antiquark production; in order to understand the origin of the different signs, however, we 
have to analyse also other diffractive final states. Contact with the corrections to the structure 



functions is made through the AGK cutting rules [14|. We find that, in contrast to the naive expec- 



tation, we cannot simply use the diffractive cross section to determine the higher-twist corrections 
to F2. There still remains an unknown piece in the initial conditions which we presently cannot 
determine. Our central result on the AGK rules is contained in eqn. ( p. 211 ). 

In the final part of this paper we present results of a numerical analysis. We write 

Fi{x,Q^) = Fr^ + ^Fi, i = t,e (1.1) 

and compare the different pieces. The twist-four correction AFi (both for i = t and i = i) consists 



of four pieces (cf. ( 2.46| )), and we will present numerical results for them, all normalized to F^ 
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Because of the uncertainty in the initial conditions, we consider three different choices, two with 
flat (in x) input distributions, and one with a rising x-distribution. As one of our main results 
for the first two scenarios, we find, despite the strong cancellations among the different twist-four 
contributions, that the total twist-four contribution is not small ai Q^ = 1 GeV^ (between 15% and 
130%). It is negative and mainly due to the transverse photon. Unfortunately, the severe limitation 
due to the DLA, together with the uncertainties of the initial conditions prevents us from obtaining 
a more precise estimate. 

The outline of this paper is as follows. In section ^ we give an overview of the higher-twist results 
which are presently available, and of the formulae that we are using in our analysis. For a systematic 
study of twist-four at small x we find that it is necessary to consider contributions of amplitudes with 
two, three and four gluons in the t-channel. We discuss these contributions in several subsections 
and outline how we handle the problem of running a^- Next (section y), we analyse the role 
of higher-twist in DIS diffraction and we study the connection with the deep inelastic structure 
function. Finally, in section ^ we describe the numerical calculations and present and discuss our 
results. 



2 Higher-Twist at small x 

We begin with a brief summary of the present theoretical status of the evolution of higher-twist 
operators. A systematic classification of twist-four operators can, in principle, be obtained using 
the arguments of R.K.Ellis et.al. 0. However, in Q (and in several other studies of higher-twist 
operators llH]) only fermionic twist-four operators have been considered. At small x, on the other 
hand, it is the pure gluonic operators that are expected to dominate, and a systematic discussion 
of twist-four gluon operators is missing. General dimensional arguments lead to the following 
operators. First, there are the four-gluon operators: 



Tr F^i° . . . F^»^ . . . F^'o, . . . F% (2.1) 



/3 

Tr i?^i° . . . p^^iP , _ _ p^^ _ _ _ pt^n 

(here I < i,j < n, and the dots denote products of covariant derivatives). In Q these operators 
have been named "quasipartonic" , and a general discussion of their evolution has been given. In 
particular, it has been shown that to leading order only two-body interaction kernels are needed. 
They are the non-forward generalizations of the AP splitting functions. An analysis of the small x 



evolution equations in DLA can be found in |16]. In addition to the quasipartonic operators ( p.l| 



there exists another gluonic twist-four operator which is contained in the BFKL equation: 

Tr F^i°D^2 . . . D^"-^F%5^^^^. + Perm . (2.2) 

In DLA the anomalous dimension of this operator can be derived from the BFKL equation; the 



recently derived |17] second order corrections to the BFKL kernel allow to go beyond DLA. In gen- 



eral, the twist-four operators (2.1) and ( |2.2| ) should mix, i.e. there should be off-diagonal elements 



in the anomalous dimension matrix. These off-diagonal elements are known in DLA-accuracy |16|. 



In the discussion below we shall see that the operators (2.1) are the most interesting ones at small x: 
they result in contributions which rise stronger at small x than the contributions of the leading- twist 
gluon operator. Therefore, at sufficiently small x, this operator will catch up with the leading-twist 
contributions, despite its 1/Q^ suppression. On the other hand, power counting in a^ will show 



that ( p.l| ) is suppressed compared to (2^): again, because of the strong small x rise in ( |2.1D , ( p.l| ) 
will eventually dominate over ( |2.2| ). The ag-suppression of ( p.ip will, however, force us to consider 
higher order corrections to ( p.2[ ). 



In the following we shall list the presently available results for the operators (2.1) and ( |2.2| ). They 
do not allow us to go beyond the DLA approximation. It remains a major theoretical effort that 
has to be made in order to calculate the full leading order evolution kernels of the gluonic operators. 



2.1 The Two-Gluon Operator 

As we have said before, the DLA approximation of the anomalous dimension of the twist-four two- 



gluon operator can be obtained from the BFKL equation. Let D^ix^Q^/Q 



t,i) denote the 



imaginary part of the scattering amplitude of the photon (virtuality Q^) off the gluon (virtuality 
Qq). It is convenient to write this amplitude as a double Mellin-transformation 



dio 



-l/2+ioo 



27ri 



du (Ql 



Dhix, Qt>/Q') = —x- :^A7^ D\{u;, u) . 



27ri VQ^ 



(2.3) 



-l/2-ioo 

In the BFKL approximation, the amplitude D\{(jj,v) can be written as 



D 



1 






(2.4) 



where D\ q denotes the quark-loop at the upper end of the gluon ladder (fig. |l|a) , and the BFKL 
characteristic function has the form 



x{^) 



Nca^ 



n 



[2^P{l)-^P{u + l)-^P{-l^)] . 



(2.5) 



It has poles at integer values of ly, and the anomalous dimensions of the leading-twist two-gluon 
operator and the twist-four two-gluon operator are obtained by solving the equation 



= w-x(i^) 
near the poles at v = —1 and v = —2, resp.: 
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(2.6) 

(2.7) 
(2.8) 



Using the recently published NLO calculations |17| it is possible to derive the singular part of the 
NLO corrections to the anomalous dimensions, i.e. corrections of the form a«^^^. 
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TABLE 1: Twist-expansion coefficients of 1)2,0 for leading- and next-to-leading twist. 



The coefficient functions are derived from the quark-loop Z)2o('^''^)- ^^ begin with the result for 
the transverse photon (^,181: 
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Similarly, for the longitudinal photon one has: 
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(2.9) 



(2.10) 



By expanding eqns. (^) and (p.lO|) in Laurent-series, 



^2,0 (^'^) 



E2 v8 Y^ 
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(y -|- n)2 z^ -|- n 



(2.11) 



and considering the residues at negative integer values of v we obtain the twist expansion in powers 
of Q\IQ'^ ■ The coefficients a„ and 6„ for leading- and next-to- leading twist are listed in table |l|. 
Since the double poles in the 6„-terms correspond to logarithms of Q^ while the single poles in the 
an-terms correspond to constants, we recognize with the help of table || that, for leading- twist, the 
longitudinal contribution is sub- leading at large Q^ (compared to the transverse cross section). For 
twist-four we have the opposite situation: the longitudinal contribution has the InQ^/Qg , whereas 
the transverse cross section is suppressed. Moreover, we note the sign structure of the coefficients 
in table y: the twist-four contribution of the longitudinal cross section comes with a negative sign. 

For our twist-four analysis we perform the computations in DLA and, therefore, consider only 
amplitudes proportional to 62 and a\ for the longitudinal and transverse contributions, respectively. 
Taking into account the evolution in Q^, we get, for each rung, a factor 72/u;(i^ -|- 2), whereF] 



72 



TT 



Performing the z^-integration we end up with 

D\{u:, Ql/Q'r=' = 4 E 4||«s (||) ' I exp (^ hqvqd) , 
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(2.12) 
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for the transverse amplitude, and 
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T7 (2-14) 



Within our conventions tlie singular part of the gluon anomalous dimension is given by 72 /tt; and not by 72 itself. 



for the longitudinal one. As stated before, the longitudinal amplitude has an extra factor IuQ'^/Qq, 
in comparison with the transverse case: in the DLA approximation, therefore, we have to be careful 
when comparing the transverse amplitude with the longitudinal one. However, the extra logarithm 
only emphasises the dominance of the longitudinal amplitude at large Q^. In the low Q^ region, 
there is no reason why the (twist-four) part of the transverse cross section should be smaller than 
the longitudinal one. We will, therefore, be slightly inconsistent with the counting of 'leading' and 
'sub- leading'. 

Finally, we note the connection between the amplitudes D2 and the twist-four corrections to the 
structure functions: 

-1/2+ioo 

-l/2-ioo 

where the superscript '2' on the l.h.s. refers to the '2-gluon' ladder. In the subsequent sections 
we will add further terms to AFi. Within our representations the twist-four part is obtained by 
evaluating the residue at z^ = — 2 (from 1/ = — 1 we get the leading-twist term). (p2{i^, ^) denotes the 
initial condition: in the context of the BFKL equation (which is an evolution equation in y = In l/x) 
1^2 would define, at fixed xq, the distribution in Q^. However, in our DLA approximation, where 
the BFKL equation coincides with the small x limit of the DGLAP evolution equation, ip2 denotes, 
at fixed scale Q^, the initial distribution in x. In the (w, i^)-representation, ip2 depends upon both 
u and LO and contains both leading-twist and higher-twist: the oo singularity structure near v = —1 
defines the x-shape of the leading-twist initial distribution of the gluon structure function, and 
the behavior near v = —2 determines the input to the twist-four correction. As usual, both input 
distributions have to be modelled and cannot be computed within pQCD. In the following, unless 
stated differently, we will use (p2 to denote the leading-twist initial distribution. 

Our discussion below will show that the situation with the two-gluon amplitude, in fact, is more 
complicated. In particular, we shall argue that the reggeization of the gluon provides extra contri- 
butions to the two-gluon operator. Namely, in the next subsection we will turn to a discussion of 



the four-gluon operators (2.1). To lowest order in as, it involves diagrams in which four gluons are 
exchanged between the quark-loop and the target proton (fig.Q). These diagrams are suppressed 
by one factor a^ compared to the two-gluon operator (^). At small x, however, this suppression 
is compensated by the stronger rise of the four-gluon operator contributions. When higher order 
gluon diagrams are included, it will turn out that, at small x, the sum of all these diagrams has to 
be split up into two separate classes. The first one corresponds to the expected four-gluon operator 
(which, however, mixes with the two-gluon operator (|2.2| )). For this piece, the ag-suppression (rela- 
tive to the twist-four two-gluon operator) is compensated by the stronger rise at small x. However, 
in addition to this contribution to the four-gluon operator there is a second set which results from 
the reggeization of the gluon. This contribution has the same form of gluon ladders as the BFKL 
equation, except that at the lower end the reggeizing gluon lines split up into two or more gluon 
propagators (fig. ^). Starting from the lower end, the four-gluon state immediately turns into a 
state of two reggeized gluons. It, therefore, has the same Q^-dependence as the BFKL amplitude, 
but counting powers of agj it is of higher order than the BFKL ladder discussed above. The lowest 
order diagram is the four-gluon diagram in fig.|2|, but in higher order it turns into a (higher order) 
contribution to the initial distribution of the two-gluon operator. We will show that this contri- 
bution can be related to the measured diffractive cross section. After the discussion in the next 

6 





z^^ 



y^ 



z^^ 




y^ 



(a) 



(b) 



FIG. 1: (a) The two-gluon amplitude and (b) D\, the part of the four-gluon amphtude which is irreducible 
with respect of reggeization. 



subsection we, therefore, will, once more, return to the two-gluon operator (^ 
separate contributions. 



and add two more 



To summarize the main results of this subsection, starting from the BFKL amplitude we have col- 



lected the known results on the twist-four two-gluon operator (2^). Rather than writing the BFKL 
amplitude as an evolution equation we have used the closed expression in the (to, z^)-representation. 
Presently, only expressions within the DLA are completely available. In order to go beyond the 
DLA it will be necessary to perform a more complete calculation. For the anomalous dimension 
in ( |2.8| ), the next-to-leading order can be extracted from the recent results of the second order 
BFKL kernel |17]. For the coefficient functions (( |2.9| ), ( |2.10| )) such a calculation has not yet been 
performed. 



2.2 The Four-Gluon Operator 

While the contributions of the two-gluon twist-four operator ( |2.2D could be derived with the help of 
the BFKL equation, the four-gluon operators ( |2.1D require a separate study. The simplest diagrams 
are shown in fig. §. In the limit of small x, these diagrams (together with higher order corrections) 
have been investigated in [^,[1^, and we simply summarize the results and relate them to the 
operators (^) and ( ^.21) . 



In the investigation of |9|,^, the starting point was the triple Regge limit of a six-point amplitude; 
from this an amplitude 1)4 was defined which describes the coupling of four t-channel gluons to two 
virtual photons. We recall only those results which are relevant for our analysis. The four-gluon 
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amplitude D4 satisfies an integral equation which we illustrate with the help of fig. |3|. As one main 
result of ref. [P, it was shown that D4 should be written as a sum of two pieces, the reggeizing terms 
Z?4', which consist of terms proportional Z?2 (eqn. ( |2.41] ) below), and terms which are irreducible 
with respect of reggeization D\: 

D4 = Df + Dl. (2.16) 

The two terms on the r.h.s. can also be distinguished by their symmetry properties: whereas D\ 
is totally symmetric under permutation of the outgoing gluons, the term D^ has mixed symmetry 
properties. In the following we shall discuss these contributions in some detail. In particular, we 
shall interpret them in terms of the operators (|2.1| ) and ( p.2| ) . In order to avoid a too complicated 
notation, we will from now on suppress the reference to the photon polarization. 

We choose our conventions such that the connection with the cross section a^*p is as follows: 



'" = -n?=w/iilJj /iiif j ^J""?". (2.17) 

In the following we will describe this formula in more detail. In particular, the overall minus sign 



in front will be derived with the help of the AGK rules |14]. 



We begin with the irreducible part. Since D^ has a particularly simple form (fig.|^ and eqn. (2.41) 
below), we simply subtract it on both sides of the equation in fig.l^ and make use of the integral 
equation for 1^2 • Regrouping the terms on the r.h.s. of fig.|^, one arrives at the following integral 
equation: 

LoDY''\kiMMM) = [D2®V'''"'^ + DY'''''®Y.^K2^2 + a-l)\ {k^MMM) ■ (2-18) 

Here fci, . . . , ^4 are the transverse components (in the Sudakov decomposition) of the momenta of 
the four gluons at the lower end in fig.|l|b, and a,b,c,d are the corresponding color labels. ^"2^2 
and a are the common BFKL kernel and the gluon trajectory function, and the sum has to be taken 
over all possible pairs of gluons. V represents the transition vertex from 2^4 gluons (fig.|l|b), 

^"'''^ " iItiI '^"''^"' [^(^l'^3)+^(^2,%) + G'(fci,fc4)+G'(fc2,fc4)- 

- G{ki, fcg + ^4) - G{k2, fcs + ^4) - G(fc3, fci + ^2) - 

- G(fc4, fci + k2) + G(fci + fe2, ks + ^4)] + 

+ 5''^5'"^ [G{ki,k2) + G{k2, feg) + G(fcl, fe4) + G{k3, fc4) " 

- G{ki,k2 + ^4) - G(fc3, fc2 + ^4) - G{k2, fei + fcs) - 

- G(fc4, fci + fcs) + G{ki + fea, fc2 + ^4)] + 

+ 5''^5'"' [G(fci, fca) + G(fc3, ^4) + G(fci, ^2) + ^(^2,^4) - 

- G(fci, fca + h) - G'(fc4, ^2 + fcs) - G(fc3, fci + ^4) - 

- G(fc2, fci + fc4) + G(fci + fe4, k2 + fea)] } • (2.19) 



The convolution D2 ® G Is given by the expression ^ 
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For the twist analysis in DLA it is useful to perform a Mellin-transformation in the following way 
(cf. eqn.(P)): 



{D2 ® G) ((J, a, b) = /-^ D2{uj, v)G{v, a, b) 



(2.21) 
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where 9 represents the angle between a and b. The leading-twist part in the DLA can be obtained 
by evaluating the residue at z^ = — 1: 
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(27r)2 ' ' 1/ + I 

whereas the twist-four contribution is related to the pole at z^ = — 2: 
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(2.24) 



Inserting the twist-two result ( 2.23| ) in ( p. 19 ), we find that the sum of the various terms cancel. 
Returning to ( p. 18]) , this cancellation means that we are loosing one InQ^, i.e. we have a non- 
leading correction to the DLA which is beyond our control and will therefore be disregarded. For 
the twist-four result, on the other hand, we find a nonzero coefficient of the pole near u = —2. 
In ( ^.18 ) the twist-four part of the ffi'st term on the r.h.s. reads (we use DLA accuracy, and we 



FIG. 2: The four-gluon amplitude at leading order. The fermion-loop symbolizes the sum over all Feynman- 
diagrams with all possible couplings of the four t-channel gluons. 



assume, for simplicity, that all kf are of the same order, k 
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(2.25) 



Here i> = 1^+2, and all momenta are normalized by Q. The factors in the first line of eqn. ( |2.25| ) arise 
from the two-gluon amplitude and the quark-loop in DLA, while the lower three lines correspond 
to the transition 2^4 gluons. A similar result holds for the longitudinal structure function (with 
a^ replaced by 621 ^iiid an extra factor 1/D in the first line). 

Making systematic use of the double logarithmic approximation, it is possible to find an explicit 
solution of eqn. ( 2.18 ). To this end we solve the integral equation by iteration: D\ then represents 
a sum of diagrams which have the following structure (fig.[l|b). At the top we start with the quark- 
loop D2,o coupled to a BFKL ladder. At the lower end of this ladder we have the 2^4 transition 
vertex V, then further below the four-gluon state with the sum over all pairwise interactions. It is. 



now, convenient to use the method of Faddeev |19]. We reorganize the sum over all pairwise gluon 
interactions and introduce auxiliary potentials T(jj)(;(.;) which describe pairs of two-gluon ladders 
between the gluon lines {ij) and (kl). The original sum over all pairwise interactions between 
the four gluon lines then translates into the iteration of these auxiliary potentials, allowing for all 
possible recombinations ("switches") {ij){kl) -^ {ik){jl) etc. For each gluon ladder there are six 
irreducible color states (8 8 = 1 © 8a © 83 ® 10 10 © 27). However, if we add the first rung 
right below the vertex in eqn. ( ^35| ) , e.g. between the gluons 1 and 2, and extract the leading 
power of hi{k'^ IQ"^) the three terms in ( 2.25| ) become proportional to ki ■ ^2^3 ■ ^4 and the anti- 
symmetric color states cancel, so that only the symmetric ones (1, 83, 27) survive. This remains 
true if one adds further rungs. As a result, our four-gluon system has only nine components: three 
color components for each coupling scheme (12)(34), (13)(24) and (14)(23). Concentrating, now, 
on the twist-four point v = —2 and collecting the maximum number of logarithms, we find the 
following pattern. For each two-ladder state described by the auxiliary potential T(y)(fc;) we have 
a propagator 7^: 






(2.26) 
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where 'KUija^ = {3, 3/2, —1} for the color states i = {1, 8s, 27}. At each switch from the two-ladder 
state {ij){kl) to the configuration (ik){jl), we collect a factor 1/2 from the angular integration, and 
insert a color recoupling matrix S. In our nine-component matrix notation, the Green's function 
S for the four-gluon system then becomes: 
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where the sub-matrices 7 and A have the following components: 
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As a result of the leading-ln Q^ approximation, at the lower end of the Green's function S the two 
ladder system Tuj\m\ ends up with the momentum configuration fcj = —kj = I and kk = —ki = I' , 
and for simplicity we set 1^ = 1' = Q\IQ'^ ■ 

For the coupling of this Green's function to the transition vertex we introduce the nine component 
vector Y\ V"^ = {v'^ ,v'^ ,v'^) with 

v^ = (2, 4\/2, 6\/3) . (2.29) 



This vector is obtained from ( 2.25D by applying the color projectors listed in the appendix. As a 
result, our function D\, written as a (transposed) nine component vector, takes the form: 
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Finally we have to couple the four-gluon system to the proton. In modelling this coupling we take, 
as a guideline, the structure of D4: to be definite, we assume that the coupling to the proton has 



the same symmetry properties as D4 in (2.16), the coupling of the four-gluon state to a perturbative 



target, e.g. to a virtual photon or to a heavy onium state. The decomposition in ( 2.16 ) reflects 
the different contributions: apart from totally symmetric (under the combined interchange of color 
and momenta) terms which are contained in D\ and -D4 , we have a term with mixed symmetry 
properties which occurs only in D^. For the symmetric piece we make the simplest ansatz, which 
contains only the color singlet representation. For the term with the mixed symmetry we simply 
follow the perturbative ansatz. This leads to: 
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(2.32) 
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and 



O fj rxz-i A/Q A/o Aj^ \ 



+ r"^ r'^'f a{1, 3; 2, 4; ^) + r^^T'^fxil, 4; 2, 3; u;) j (2.33) 

with positive- valued functions /s and /a- The signs follow from the AGK rules and will be discussed 
below in more detail. The pre-factors 1/8 and 1/3 are convenient for color and for statistics, resp. 
We immediately see that the (/94A-term does not couple to D\ but only to D^, and we will come 
back to this term when we discuss -D4 . In order to discuss the momentum structure of the ansatz 
(/74s, let us recall that in the usual DIS leading- twist ladder the lowest cell has only a longitudinal 
integration, and the integration over the virtuality of the lowest parton is absorbed into the initial 
condition lp2{Q'q,x). In evaluating our four-gluon system in the DLA we treat each gluon ladder 
of our auxiliary potential Tuj\m\ in the same way as the leading- twist ladder. Let the lowest 
potential be 7'(i2)(34)- At its lower end we have the momentum factor [kik2){k^k4). In the function 
fs{iij] k, I; io) of our ansatz ( p. 31 ) the main contribution comes from the region where 



ki = I + r, kj = —I + r, k^ = m — r, ki = —m — r, 

r^ <c /^ m^ . (2.34) 

Combining these functions with the momentum factors from the auxiliary potential T'(x2)(34) we are 
lead to the following definition of our initial condition: 



o f'^^dP f^odm 



,2 



Qo¥'4s(Qo;^) = / dr / -wl — ^l m fs{ki,kj;kk,ki;uj). (2.35) 

Inserting into (|2.32| ) the momentum assignment ( p.34| ), and using the definition ( |2.35| ) we get a 
factor 

where the factors 1/2 in front of the second and the third term are due to the angular integrations 
over I and m (a similar definition holds for the functions 994A and /a). As to the longitudinal 
integrals, since we work in w-space we allow the function (/94s to depend on a;: if we assume a single 
ladder to rise as (1/x) , 9343 should have a pole at uj = 2A, i.e. the initial distribution rises twice 
as strong as the leading-twist gluon structure function at the same scale. 

Finally we couple this ansatz to D\ in ( 2.30 ). Projecting onto the color eigenstates and introducing 
the nine-component vector 



*4s(Qo,'^)=h^ •'^4s(Qo;'^), "" = ^7^272 (2.37) 

we arrive at the following expression for our twist-four contribution to the structure function in the 
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FIG. 3: The integral equation for D4. The sums denote couphngs to the gluon hues in all possible ways. 
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and an analogous expression in the longitudinal case. After performing the inversion within S = 
{G~^ — ^S)~^ we get explicit expressions for the components of vector T,V. They are listed in the 
appendix. 



As one of the most striking features of ( 2.38| ) we note that S has a leading singularity in the P-plane 

at m,!^ 



D = 4{l + 5) 



iVrOs 



TTLO 



(2.39) 



where 5 = 0.009549. Using the usual saddle point approximation, this pole leads to the following 
asymptotic small x behavior: 



^4 ~ ( ^)'exp (2 (1 + ^^aMiMHQVQl)] ■ 



(2.40) 



This rise at small x is approximately twice as strong as that of the leading-twist gluon structure 
function, and at sufficiently small x the four-gluon higher-twist term will become as strong as the 
leading-twist term, despite its suppression coming from the extra power in ag and the 1/Q^ factor. 



This strong rise at small x is the reason why the four-gluon operator, (2.1), is expected to be so 
important in the small x, low Q^-region. 



Besides the leading pole ( p.39|) , there are more singularities in the P-plane, which have been dis- 
cussed in some detail in ref. [M|. There are three cuts in the P-plane. The first one stems from the 



singlet-part of matrix 7 (eqn. (2.2S|)). It starts at t' = and goes up to the point u = 4 ^"'^ which 



is very close to the leading pole ( |2.39| ). It overlaps with the symmetric color octet cut, which is 
located in the interval [0, ^]. Finally, the cut of the color 27plet goes from — ^ to 0. 
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FIG. 4: The reggeizing part of the four-gluon amphtude has the same evolution as the two-gluon ampUtude 
but is suppressed by one factor Us- 



Finally, we return to (|2.16| ) and we discuss the reducible part D^. From |9|,|l6| we have: 

Df{ki,k2,k3M)= (2.41) 

5 I ^afccd^^^^^^^^^^^^^^^)^^^^^^^^^^^^^^^^ _ D2{ki + kA,k2 + k3)] + 



2^2 



+ d'''"^[D2ik2, fcl + fcs + ^4) + 1^2(^3, fcl + fc2 + ^4) - ^'2(^1 + ^2, fcg + ^4) 



D2{ki+k3,k2 + k4)] }, 



where the color tensor d"'^'^'^ is given in the appendix. It contains both symmetric pieces and 
pieces with mixed symmetry: one easily recognizes the structures ( |2.32| ) and (2.33). Convoluting 
D^ with our ansatz ( 2.3l| ) we find that ip^s couples only to the symmetric and ip^A only to the 
mixed-symmetric part of D^. With the help of the same arguments as before, cf. ( 2.35 ), we arrive 
at 
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(2.42) 



Here (/34A is defined in analogy with (/94s in ( 2.35| ). It is important to note that the negative sign 
in front of (/?4a is dictated by the AGK rules (see below). Note that, to lowest order in Og, it 
is this contribution u^ that contains the diagrams with four elementary gluon lines between the 
quark-loop and the proton (fig.^). 

Having written down our DLA expressions for the four-gluon twist-four contributions, we briefly 
return to the operator product expansion. As before, the ladder diagrams denote the In Q^ evolution 
in the DLA, but instead of writing the four-gluon amplitude as an evolution equation in InQ^ we 
have used the closed expression in the (w, i/)-representation. Twist-four corresponds to the point 
V = —2. D\ then has the following interpretation: the upper two-gluon ladder (evaluated near 
u = —2) describes the evolution of operator ( |2.2| ), the 2 — > 4 vertex the mixing between ( |2.1[ ) 
and ( |2.2[ ), and the Green's function S the evolution of the four-gluon operators (2.1). Obviously, 
the three auxiliary potentials T'(i2)(34), 7'(i3)(24)) ^{14)(23) belong to the three operators in the first, 
second, and third line of (p.lj), and what we have called "switching from {ij){kl) to {ik){jl)" has 
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to be interpreted as the mixing between these three different four-gluon operators. If we were to go 
beyond DLA, there would also be a transition from the four-gluon state back to two gluons (this 
vertex could again be derived from ( 2.19| ) and ( 2.22D , by expanding about z/ = 1), and eventually 



these transitions would be iterated and exponentiated in the i-channel. A remarkable feature of 
the small x approximation (upon which the DLA is based) is the fact that the four-gluon operator 
has no direct coupling to the quark-loop. Our ansatz ( p.31| ) for the coupling to the proton, i.e. for 
the initial condition, which we so far have motivated by inspecting the structure of -D4, also arises 
from an analysis in the spirit of Q : the matrix element of the product of four-gluon fields has to 
be decomposed into terms with different symmetry properties. 

The interpretation of the contribution D^, on the other hand, is more subtle. As stated before, this 
term has its origin in the reggeization of the gluon: from the point of view of DIS the appearance of 
reggeization is a "novel" phenomenon since at the leading-twist level it has played only an indirect 
role (e.g. as part of the NLO anomalous dimension of the gluon). From the point of view of Q^ 
evolution, D^ belongs to the two-gluon operator (p.2|). Including D^, therefore, at first sight 
seems like adding a non-leading correction to ( p. IS] ), and within the DLA we better should restrict 
ourselves to the leading term (|2.15| ). On the other hand, the analysis in [^ clearly shows that D\ 
and D^ have the same origin (diagrams with four gluons in the t-channel) and come with the same 
powers of logarithms. We, therefore, believe that for the twist-four calculations - which necessarily 
have to include t-channel states with four gluons - we have to keep D^ as a separate contribution. 



in addition to ( 2.15 ). 



A better understanding of why we have to keep this contribution can be obtained by looking at 
the evolution of the gluonic operator in ( |2.2[) . The mixing between ( p.l| ) and ( p.2| ) is of order a1 
in the anomalous dimension matrix, i.e. in total the correction to the evolution of (|2.2| ) due to the 
transition from ( |2.2| ) to (2.1) and back to ( |2.2| ) is a NNLO effect of order a^. To be consistent, 
we, therefore, have to include NNLO corrections also in the anomalous dimension 7^^"^ in 



and in the coefficient function. They are not available yet. On the other hand, what we are seeing 
in D4 can be viewed as a particular piece of these corrections, connected with the reggeization of 
the gluon. Moreover, there are reasons to expect that they are the most important ones. Since 
we are studying t-channel states with four gluons, D^ provides contributions from the "decay of a 
reggeizing gluon into two or three gluons" which then become partons inside the proton. In other 
words, we are encountering the two- and three-particle states in the gluon trajectory function, 
appearing at the lower end of the two-gluon ladder. Clearly, these higher-order particle states in 
the trajectory function can appear not only at the lower end of the two-gluon ladder (as taken 
into account by D^) but also somewhere between the proton and the photon. Such a contribution 
would then be part of the NNLO corrections to the gluon rung, i.e. to the anomalous dimension 
7"^^^ in (|2.S| ) mentioned above. In this sense, D^ counts those NNLO corrections to the two-gluon 
kernel which are connected with four-gluon states in the t-channel, and these corrections are taken 
into account only at the lower end of the ladder. When running Og and the strong ordering of 
the transverse momenta is taken into account, we expect these corrections to be largest when they 
appear at the lower end of the ladder. For these reasons our incomplete treatment of the NNLO 
corrections in the evolution of the two-gluon operator ( |2.2| ) may not be such a bad approximation. 
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2.3 Contributions with Three t-channel Gluons 

Since there are no three-gluon operators that might contribute to our twist-four analysis of the 
unpolarized structure functions, t-channel states with three gluons can come in only through the 
mechanism which we have discussed at the end of the previous section. In the evolution equations 
of the two-gluon operator ( |2.2| ) we should include the NLO corrections to the anomalous dimension 
7'^^'' in (|2.8D as well as to the coefficient function. Whereas the former ones are, at least in principle, 



now available (so far, they have not been computed from [17|; as pointed out in |2^, it is not clear 
how the behavior near u = —1 can directly be extracted from the results of [0), NLO corrections 
to the coefficient function are not known. In our analysis, therefore, we follow the same logic as 
for the four-gluon case. Prom ^ we know scattering amplitudes with three gluons in the t-channel 
(in in they are called -D3). Because of the reggeization of the gluon they can also be written as a 
sum of D2 functions (cf. ( |2.4lD ): 

D3(fcl, fc2, fes) = ^r'' (^2(fcl + fc2, fcs) + D2{ki,k2 + fcg) " 1^2(^1 + ^3, ^2)) • (2.43) 

The convolution with the initial distribution of the proton requires a new function, c/?3(w). One ends 
up with functions AF.^ ^-, that have the same dependence on variables as ( |2.42| ) (with a different 
constant factor in front of the r.h.s.; it will be absorbed into the unknown initial distribution 993). 

2.4 Running as 

Since in DLA the leading-twist two-gluon amplitude must coincide with the DGLAP solution, it is 
clear how one can consider the fact that Os is running in this case. One has to exchange the value of 
fixed as by Og = 47r//3o (72 — > 12//3o) and the expression Q'^/Qq by t/to, where t = ln((5^/A^) and 
to = ln((5o/A^). This procedure should be applicable for the bare four-gluon amplitude (i.e. the 
amplitude without two-gluon amplitude and vertex at the top of fig.|l]b), as well, since the addition 
of one rung is always performed by convolution with the common BFKL kernel. This means that 
every time as comes with a factor In (t/to) it has to be replaced by ag. 

There are three powers of ctg in eqn. ( |2.38| ): One stems from the quark- loop, which (in the longitu- 
dinal case at leading order) is proportional to a logarithm. The remaining two powers stem from 
the coupling of the four-gluon to the two-gluon amplitude via the Vertex V (fig. ||b) . Here we get 
only one logarithm. We write this coupling symbolically as 

dY = £dr'exp[72(r - T')]fa', |^ expi^r' - ro)]G(^) , (2.44) 

where r = Int etc., G(z^) is understood as the four-gluon amplitude Greens-function and / shall 
contain all remaining factors. In eqn. ( 2.44| ) one factor Og needs to be replaced by ag(t') and the 



other one by Qg. If, on the other hand, we go back to fixed Og, we have to replace r and r' in 
eqn. ( |2.44D by t and t', resp. Performing the t'-integration and comparing with eqn. ( |2.38| ), we can 
identify the quantities G and /. By this way we arrive at the following conclusions: Por running Os 
the factor {Q"^ /Q'^Y ™ ^qn. ( 2.38| ) has to be replaced by {t/t^Y , one power Og has to be evaluated 
as as((5^), while (in the longitudinal case) the remaining two powers are replaced by Og. In the 
denominators we have to replace the factor {uju — 72) by {uju — oo — 72). 
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FIG. 5: Two of the four diagrams which enter into diffractive qg-production at leading order. 




FIG. 6: The twist-four part of diffractive qg-production is described by exchange of gluon ladders. 



Repeating in ( ^.301) the standard saddle point analysis (cf. eqn. ( ^.391) ), we arrive at a small 
behavior of the form 
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(2.45) 



Apart from the small 5-correction, this small x behavior is just the square of the well-known leading- 
twist double scaling formula ||23|,p^]. Therefore, in the kinematical regime in which D\ gives the 
main contribution to twist-four, higher-twist is expected to increase much faster with decreasing x 
than leading-twist and to decrease slower than 1/Q^. 

We end this section with a brief summary. All in all we have collected four different contributions 
(both for the transverse and the longitudinal structure function). In our notation 

AF = AF(2) _^ ^^(3) _^ ^^R _^ ^^i ^ ^2.46) 

where we have suppressed labels t and i. As far as the Q^-evolution is concerned, the first three 
terms belong to the twist-four two-gluon operator (|2.2|) whereas the last one belongs to the four- 



gluon operators (2.1). This last contribution is of particular interest, since at small x it rises stronger 
than the other ones and, therefore, may potentially become large. However, we have already seen 



that the four contributions in ( 2.46 ) come with alternating signs, which hints at the possibility of 
strong cancellations. 



3 Diffractive Dissociation 

Up to now, we have been describing the four contributions to twist-four, which we believe to 
dominate at small x: D^, D\, D^ and the twist-four part of the two-gluon amplitude. Next, we 
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Contribution 


Sign 


HQVQl) 


/3^1 


Trans, r = 2 


+ 


non-pert. 


vanishes 


Trans, r = 4 


— 


no 


vanishes 


Long. T = 2 


vanishes 


vanishes 


vanishes 


Long. T = 4 


+ 


yes 


constant 



TABLE 2: The four contributions to difFractive gg-production up to twist t ~ A enter with different signs into 
the diffractive cross section. The table shows also which contributions have leading logarithms, 
and how they behave for /? — > 1. 



have to address the important question of specifying the initial distributions tp. Following the 
conventional leading-twist analysis, we would determine the initial conditions from a fit to the 
data. In this first attempt of analysing the influence of twist-four corrections at low Q^ and small 
X, we will try to estimate the ip functions with the help of empirical data on final states which are 
known to belong to twist-four. The best candidates for this are diffractive final states, in particular 
the production of longitudinal vector particles and the production of jets with large transverse 
momenta. In this section we, therefore, give a brief overview over the role of higher-twist in DIS 
diffraction. In particular, we will discuss the relation between the higher-twist pieces of diffractive 
cross sections and the higher-twist contributions to Ft and Fi from the previous section. In our 
numerical analysis these relationships will be used to estimate the initial distribution (/'4S,A- 



3.1 Diffractive qg-Production 



It is useful to divide the kinematical regime into different regions of the variable j3 = Q'^ / {M'^ + Q'^) , 
where M is the invariant mass of the diffractive system. We begin with the small mass region 
(/? — > 1), where the gg-production process is expected to give the main contribution to diffraction. 
Those contributions to the cross section, which can be calculated perturbatively, can be described 
as follows. Seen from the angle of the proton rest frame the virtual photon splits up into a qq-pair, 
which interacts with the proton before the quarks result into the two jets. For diffraction the 
interaction between the qq-pair and the proton must be colorless and, therefore, in the small x 
regime where gluons are expected to dominate, at least two gluons need to be exchanged (fig.|5[). 
The cross section can, now, be computed by coupling the gluons in fig.|| to the proton and by 
squaring the diagrams afterwards. This leads to the following expressions in the case of zero 
momentum transfer (t = 0) [p5, 26|: 
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Here pt is the transverse part (in the Sudakov decomposition) of the quark-momenta, Xjp = x-q/ 13-, 



2 ;2\ 



l\Q\M\pi,l 



2 J2 i2\ 



r{Q\M\pi,i 




and Fq represents the unintegrated gluon distribution of the proton 



(3.3) 



(3.4) 



(3.5) 



In order to estabhsh a connection with the deep inelastic structure functions we have to integrate 
over M^ (or (3), pf and t. Performing in eqns. ( |3.1| ) and ( |3.2D the pf- and M^-integrations one 
arrives at the quark- loop at the top of fig. g , with the constraint that the two-gluon systems at the 
r.h.s. and at the l.h.s. are in a color singlet state. We denote this with the help of the superscript 

(l,+,+):rfo^'+). 



In terms oi D\ n ' the integrated diffractive cross section becomes: 
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(we have suppressed the distinction between longitudinal and transverse polarization). As shown in 
P, D\ q ' (fci, k2, ks, k^) can be written as a sum of D2,o expressions with appropriate evaluation 
of momenta: 
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4,0 



'^\ki,k2,k3,ki) = g'^ — { 1)2,0(^1,^2 + ^3 + ^4) +-^2,0(^2,^1 + fcs + ^4) + 

+ ^2,0(^3, fci + ^2 + k4) + D2fi{k4, ki + k2 + ks) - 
- D2fl{ki + k2,k3 + ki) - D2fi{ki -F fes, ^2 + ^4) - 
-D2fl{ki + k4,k2 + k3) } . 



(3.7) 



Here fei , . . . , ^4 are the transverse components of the gluon momenta at the lower end of the quark- 
loop. In our case the quark-loop is coupled to the proton through the unintegrated gluon structure 

V 



function !Fg (cf. ( |3.lD and (3^)), and we have to substitute fci = — ^2 = I and ^3 = —k^ 
and integrate over I and I' . To define leading- twist and twist-four corrections we simply expand in 
powers of 1/Q'^. 

Before we describe the formal expansion in inverse powers of Q^, let us give a qualitative descrip- 



tion |27|. We begin with the p^-integral of the transverse cross section ( ^.1[ ) and keep M^ fixed. 
As long as p^ is not small, the leading contribution of the Z-integral comes from the region of 
small P < p1{M'^ + Q'^)/M'^ where T^ /P behaves as a constant. The Z-integral, therefore, simply 
leads to x^g{x^,p1{Q'^ -j- M'^)/M'^), and the transverse cross section falls as [l/pt x^g{x^,p'^{Q'^ + 
M'^)/M'^Y- Using naive dimensional arguments it follows that this region belongs to higher-twist. 
Taking the integral over p^ one finds dominance of the low pi region where perturbation theory 
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breaks down. The leading-twist of the transverse cross section is obtained by extending the conver- 
gent p^-integral up to infinity. The twist-four term arises as a correction: the correct upper hmit 
to the p^-integral is Q^(l — /3)//3, i.e. we have to subtract the integral from Q^{1 — 13)/ 13 to infinity. 
This is the negative twist-four correction to the transverse cross section. The convergence of the 
Pj -integration implies also that there is no In Q^ in the twist-four correction to the transverse cross 
section. Turning to the longitudinal cross section, one notices in ( p. 21 ) the extra pi factor in front 
of the square brackets: this changes, compared to the transverse case, the p|-behavior of the cross 
section in two ways. When integrating over pi, the small p^ region no longer dominates, whereas 
in the large p|-region we encounter a logarithmic divergence. This explains the 1/Q^-suppression 
of the cross section and the appearance of a InQ^. The (3 (or M^) dependence of the two cross 



sections arises from a closer inspection of the two formulae ( |3.1| ) and (3^): whereas the transverse 
cross section (both leading- and higher-twist) vanishes near /? = 1, the longitudinal one stays finite. 
This leads to the conclusion that the longitudinal cross section (which for /3 <C 1 is much smaller 
than the transverse one) dominates in the large /3-region. In other words, the diffractive qq cross 
section near /3 = 1 is mainly longitudinal and belongs to twist-four. Results of this discussion are 
summarized in table ^. 

In order to obtain quantitative expressions for the higher-twist cross sections (and to verify our 
intuitive arguments), we return to the p^ and /3- integrated expressions in (|3.7| ) and perform a formal 
expansion in powers of 1/Q^. Using the (a;, i^)-representation for 1)2,0 iii ( P-'^ ) and expanding about 
the twist-four point u = —2 yields for the transverse case 



<+'+)(., i\eY=^ = 4^i Y. 4^«s {2/^ + 21" - {I + ir - {I - ir} ■ m 



f 

Similarly, in the longitudinal case we get 



il)lJ+'+)(^, I', fY=^ = -hi^l E 4^°'-^"" '" " ^ "'" ^" '" 



/ 



{I + I'f \n{l + I'f -{I- I'f ln(Z - I'f^ . (3.9) 



Here, /^ and /' are normalized by Q^. Working in DLA, we have for the unintegrated gluon 
structure function !Fq, an expression of the form 



,2, 1 d fduj fl\ /72i„,2/ 2\ ^ /, ,Nr=2 



^o(-. I') = Qjjp y 5- [-) exp [^ 1„ l'/,i) :M-Y-' . (3.10) 

where '^2{<^) is the initial distribution of the gluon structure function and q^ = Qf^/Q'^- Multiplying 
with functions ^(^{1"^) and J-qQ.' ) and performing the integrations over I and V , we arrive at 
the following higher-twist correction to the differential diffractive cross section: 



A 



dt 



t=o 



^^5^e^exp(2^1nQVQ^) [MQI-^)? ■ (3.11) 



Q' 6 .. ^ 



For the longitudinal case the Z, V integrals are slightly more complicated, since ( |3.9| ) contains 
logarithms. In order to obtain the maximum number of logarithms, we split the integration interval 
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FIG. 7: Duality relates the diffractive vector meson production cross section to the difFractive qq-production 
cross section. 




FIG. 8: The (integrated) ggg-production amplitude is given by an expression which is symbolized by this 
figure. Note that the effective 2-^4 transition vertex is not the same as in fig.|l|b for D\. 



into subintervals /^ < /' and /' < l^. The result is 



A 



dfjOT"' 



dt 



47r^ap 



at 



t=o 



Q4 67r2 272 



Y,e}{l-exp{2^lnQyQl)][MQl^)?- (3-12) 



Before we can draw the relation with the inclusive structure functions we still have to integrate 
over the momentum transfer t. Repeating the discussion in the sequel of eqn. ( |2.33| ), we note that 
in DLA the nonzero momentum transfer enters only into the lowest part of the ladders, the initial 
distributions ip2- Integration in ( |3.11| ) or ( p. 12 ) over t, therefore, simply means replacing (/j^ by 
f d'^q(p2, i.e. by an effective new initial condition. Moreover, if we allow for diffractive dissociation 
of the target proton, this initial condition will again be modified into a new initial condition which 
we denote by 994'^ (w): 



Aa'?9/ 



An'^ap 



at 



Q2 67r2 272 






Ql 



exp 



2^1nQVQo 






(Ql^ 



(3.13) 



(and a similar expression for the transverse cross section). Further below we shall discuss how this 
diffractive cross section contributes to the structure functions Fi^f In particular, we will derive a 



relation between ip'^^^{uj) and (pis,A{i^) introduced in ( 2.3l| ). 
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FIG. 9: Same as in fig-g] but now for qqg g-production. The sum denotes all possible couplings between 
systems (12) and (34). 



3.2 Diffractive Vector Meson Production 



Let us, now, turn to those diffractive processes whicli have been examined most intensively both 
theoretically and experimentally, the diffractive production of vector mesons: 7* + p — > y + p, 
where V can be any vector meson. For the longitudinal photon this process has been shown to be 
calculable within perturbative QCD |10|-13|: 



da. 



'y*p-^Vp 



dt 



t=o 



12TT-^ai My 



N^ae 



V- 



-TiQ'H 






d 



2 dlnx 



xg{x,Q'^ 



(3.14) 



Experiments show that the sum over all vector particles gives a significant (about 20 %) contribution 
to the total diffractive cross section. As can be seen in ( p.l4| ), for an individual vector particle, 
the cross section goes as 1/Q^. In ref. [13| it has been shown that the production cross section 
of a single vector particle can also be obtained if one starts with open gg-production and projects 
onto the corresponding angular momentum and parity quantum numbers. This suggests to apply 
a simple duality argument and to set, for the mass range rup < M < Mc, the sum over vector 
particle production cross sections ( |3.14 ) equal to the M^- (and p^-) integrated cross section of 
open Q'^-production ( |3.2[ ) (we illustrate this equality in fig.|^. In this way, the sum over the vector 
particle cross sections turns into a twist-four contribution: the additional factor l/Q^ in ( 3.14 ) is 
due to the projection onto the vector particle wave function. As a result of this argument, we can 
use the sum of the measured vector particle cross section to obtain a lower bound of the integrated 



longitudinal qq cross section |28|, which has been shown to belong to twist-four. 



3.3 Diffractive gg(7-Production and Multi-Jet Processes 

Diffractive (/^-production is expected to be dominant near /? = 1. If we go to /?< 1, we have to 
consider also gluon production processes {qqg-, qqgg-,- ■ .-jets). Cross sections for these processes 
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Contribution 


Sign 


HQVQl) 




Trans, r = 2 


+ 


HQVQl) 


non-pert. 


Trans, r = 4 


+ 


HQVQl) 




Long. T = 2 


+ 


no 


non-pert. 


Long. T = 4 


- 


{HQ^/Ql)f 





TABLE 3: The four contributions to diffractive qq(7-production. We list the sign structure and powers of 
lnQ2. 



have been calculated in different kinematical regions: for very small (3 (the triple Regge region) 
in 1^,29], for strong ordering in the transverse momenta in |3C, 31|. The latter calculation allows 



to describe the whole /3-interval, but only the first one contains both leading-twist and twist- 
four. In this section, therefore, we make use only of the results obtained in |2£] and discuss the 



higher-twist contributions to diffractive g^^f-production. The calculations are done in analogy to 
the gg-production process discussed above: we keep all powers in the transverse momenta, but we 
restrict ourselves to the leading log in M^ (or 1//3). An example of the Feynam diagrams is shown 
in fig. |To|a. A complete description can be found in [^. In general, the s-channel gluon can be 
emitted anywhere from one of the t-channel gluons. 



The analytic expressions for the cross section formulae [^9| are to lengthy to be repeated here. We 
follow the same logic as in the sequel of (^]^) and (|3.2| ) . After integration over the momenta of the 
outgoing quarks and the gluon the result |^ takes the form illustrated in fig.|^: all diagrams can 
be drawn in such a way, that at the upper end we have the fermion-loop -D2,0) ^-nd the s-channel 
gluon is contained in an effective vertex function which is similar (but not identical) to the 2 — > 4 
gluon vertex in ( ^.191) (fig.||b). From p we take the expression: 

^mflil^ I') „ = 14^2,0 ® {2G(i, -I) + 2G{1', -I') + 2G{1, 1') + 2G{1, -I') + 



t=o 4\/2 



+ G{1 + I', -I - I') + G{1 - I', -I + I') - 2G{1, -I + I') 

- 2G{1, -I - I') - 2G{1', I -I')- 2G{1', -I -/')}. (3.15) 



Here G{a,b) is the same function that enters into the convolution D2 V (eqns. ( p.lSj) , ( |2.19| )). 
The ladders at the lower end of this vertex, again, denote (leading-twist) gluon structure functions 



Since we did not present the explicit formulae for the (unintegrated) diffractive cross sections, we 
only briefly sketch the qualitative picture (details can be found in [p9l). Let k2 be the transverse 
momentum of the gluon, and fci, — (fci + ^2) those of the quarks. The easiest way is to start 
with the hard region (i.e. the transverse momenta of quarks and gluons are large), and to consider 
configurations where the momenta are ordered: fc^ < kf. In this region, the transverse cross 
sections goes as da* ~ 1/[A;^A;2] • [Xpg{x^, k"^)]"^ and the longitudinal cross section behaves as dcr^ ~ 
l/i^i^f] ■ [^]p5(^]p) ^i)]^- ^y dimensional counting, this region belongs to twist-four. Leading-twist 
can be obtained if one realizes that the integration over the gluon momentum ^2 is dominated by 
the small momentum region, i.e. by nonperturbative physics. The remaining fci-integral provides, 
for the transverse photon, a InQ^-enhancement, whereas for the longitudinal photon it does not 
(table I). 

The derivation of the twist-four contributions, now, follows the same way as for gg-production. 
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We use the (a;, z^)-representation for ( |3.15| ), and the twist-four part (near v = —2) takes, for the 
transverse photon, the form: 






4i ,/2 






TT W 



2^4 + 21'^ - 8^2/'^] In [max(Z2, Z'^)] } . (3.16) 



Similarly, for the longitudinal photon we obtain: 



2 7/2N,r=4 



W 



^L-,,o('^,^^^') 



b\ 9 1 y- 



2 vr a; 






,2 ^,3 Q 



f 



Qo •- 

21^ + 2/'^ - ^IH'^] In^ [max(Z2, r^)] } . (3.17) 



The integration over I, V is done in the same way as described after ( |3.10| ). With ( p. 10 ) for the 
unintegrated structure function we find 



A 



do-'J^s-t 



dt 



47r^aem San 



t=o 



Q4 167r3 27; 



272 a; ^ ^ [ 



1 — exp 



272, Q' 



^l{u) 



(3.18) 



and 



A 



do-959.^ 



di 



47r^ap 



t=o 



Q4 32 



^3 2^2^ /\ 



exp 



272, Q^ 

u Qf. 



272, Q^ 



(/9^(tj) 



(3.19) 



for the transverse and for the longitudinal photon, resp. Finally, performing the integration over 
t and replacing ip^ by Q'^'^'^^ leads to our final expressions for the twist-four corrections to the 
diffractive structure functions Ao"'^'^^'* and Aa'^^^'^. 

An important feature of these results is the sign structure (summarized in table |3|). For both 
the transverse and the longitudinal photon we find the opposite signs compared to gg-production 
(table pi). This implies that when we are calculating the total twist-four corrections to the diffractive 
cross section - either transverse or longitudinal - we are adding two terms (gg and gg^f) of opposite 
signs. Moreover, the higher-twist corrections for the transverse and for the longitudinal diffractive 
cross sections have opposite signs. This clearly allows for the possibility of substantial cancellations, 
in particular in the corrections to F^ = Fp + Fp. The longitudinal diffractive cross section may, 
therefore, be a cleaner place to look for twist-four corrections. 

It is not difficult to generalize this discussion to the production of more gluons in the region of 
very small /3. For example, for the (integrated) production cross section of two gluons we have to 
calculate the diagrams shown in fig.|^ As before, use has been made of eqns. ( p. 41 ) and ( |2.43| ). 
In this way all diagrams in fig. ^ can be rearranged in terms of the vertex ( 3.15| ), and they can be 
grouped into two classes. Beginning with the fermion-loop at the top of the diagram, we have, in 
the case of fig. ^, one BFKL rung, then the vertex ( |3.15| ) , and further below the two-gluon ladders 
for the unintegrated gluon structure functions. In case of fig.^a, there is no rung between the 
fermion-loop and the vertex ( |3.15 ), but instead four possible rungs below the vertex at the top of 
the gluon ladders. Similarly, for three gluons we have three possibilities: two gluon rungs above 
the vertex, one rung above and one set of rungs below, or two sets of rungs below. 
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3.4 Connection with the Structure Functions F^i and the AGK Cutting Rules 

After we have discussed the different diffractive cross sections, we, now, turn to the question of 
how these contributions enter into the inclusive cross section ( p.l7 ). An important ingredient are 
the AGK cutting rules ||l4|. 



We first return to the total cross section ( |2.17 ) and recapitulate how the AGK rules work. Let us. 



for the moment, replace the proton by a virtual photon, i.e. we consider the elastic scattering of a 
virtual photon on another virtual photon. This process can be treated perturbatively. For the time 
being we do not expand in powers of Q'^/Q'^-, but consider the Regge limit at large but finite Q^. 
In order to study the AGK rules, we start with the discontinuity across the four-gluon t-channel 
intermediate state. Disregarding all unnecessary details, this discontinuity takes the form 

disc^ F^.^. ~ -CI C4 (3.20) 

where F^*^* denotes the t-channel partial wave for the 7*7* scattering process, C4(A;i, A;2) ^3, ^4; ^) = 
{uj — ^ (i{ki))Di{ki,k2-,k2iki'-,oj) is the amputated 7*7* — >4-gluons partial wave, and the symbol 
contains all phase space factors. For C4 we have the same decomposition (eqn. (|2.16)) as for D4,. 



In particular, C\ is completely symmetric, and C4 contains a symmetric piece and a piece with 



mixed symmetry (cf. the discussion after (2.41)). In the unitarity integral (3.20), we take from 
both C4 and C| either the symmetric pieces or the pieces with mixed symmetry, i.e. there is no 
interference term. 



First consider the symmetric pieces, e.g. C\ for both factors in ( 3.20| ). Writing C\ in terms of 



the auxiliary potentials T/ij\ii^n from our discussion above, we have in our unitarity integral on 
the r.h.s. of ( |3.20|) the three identical 'diagonal' terms T'(;^2)(34) ^ 2^(i2)(34)i ^(13)(24) ^ 2^(i3)(24)i ^^'^ 
^(14)(23) ®2^(i4)(23)- Together with the statistical factor 1/4!, this gives a weight factor 1/2(1/2)^ for 
the diagonal term (this counting does not include color). Similarly, 'nondiagonal' terms, Tn2)(34) ® 
^(13)(24) 6tc., yield the weight factor (1/2)^. In order to verify the AGK cutting rules we compare 
2Imr^.^* of this amplitude with the different cuts, (Tq, cti, and (T2 (the subscripts refer to the 
number of cut ladders) . The (diagonal) "diffractive cut" uo has the cutting line between the gluon 
lines '2' and '3'. Its contribution is positive and has the weight (1/2)^ (due to the statistical factors 
inside the ladders). For the "double multiperipheral" cut <T2 the cutting line runs also between 
gluons '2' and '3', but its (positive) contribution has the weight 2(1/2)^. Finally, the "absorptive 
cut" o"! is negative; the cutting line runs between line '1' and '2' or '3' and '4', and the counting 
gives —2 • 3 • 1/3! = —1. Adding all these contributions, we arrive at —(1/2)^, i.e. minus the 
diffractive cross section ai. The result agrees with 2ImT^.-y*, as anticipated by the AGK rules. 
Two features of these results are important for us: in total, the four-gluon t-channel state gives a 
negative correction to the single-ladder contribution (this justifies the minus sign in ( 2.17| )), and 



its absolute value equals the 'diffractive cut' ctq- It is easy to verify that the same conclusion also 
holds for nondiagonal pieces T'(]^2)(34) "25 Ti^i3)(24) ^tc. 



Next, the part with mixed symmetry: in ( 3.20 ) we put Cf' for both C4-factors, and we take the 
terms with the anti-symmetric structure constants (|y™c jccd Q^Q^y Their color and momentum 
structure is the same as in ( |2.33D . They are connected with the odd signature reggeized gluon. 
Again, in the unitarity integral we have diagonal terms of the type /a(12,34) ® /a(12,34). Their 
total weight is 2(1/2)^. In the counting of the nondiagonal terms we encounter several cancellations: 
the weight factor turns out to be 1/2. Now, we compare this with the cuts uo, <ti, and (T2. For the 
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diagonal (Tq and (T2 contributions where the cutting hue runs between gluon '2' and '3', we find the 
weight factors (1/2)^ and 2(l/2)'^, resp. They are the same as the counting factors for the even 
signature ladders which we have discussed before. In particular, we again have the negative sign 
for the sum of the three cuts. There is, however, an important difference between even and odd 
signature reggeons: for odd signature the real part of the signature dominates, and the leading 
contribution cjo is contained in the LO BFKL equation. The only contribution to o"o from the four- 
gluon intermediate state comes from the 'diagonal' pieces which present the two gluon contributions 
to the gluon trajectory functions. Consistency, therefore, requires that the magnitude of this term 
is linked to the leading order BFKL ladder. One can check from ( ^.4lD that this condition is, 
indeed, satisfied. 

In summary, the AGK rules allow us to determine the sign of our initial conditions. We return 
to deep inelastic scattering off the proton, consider the twist-four term in the expansion in powers 
of Qq/Q^ and replace, in ( 3.20 ), the second C4 factor by our general ansatz ( 2.31 ), i.e. by our 
nonperturbative initial conditions V34s and ip^A- In order to preserve the sign structure dictated 
by the AGK rules, both functions have to be positive. Moreover, since (/94A is connected with 
the reggeization of the gluon trajectory function, it is linked to the size of the twist-four part in 
the BFKL ladder, /S.F^'^K In this first attempt to estimate the higher- twist contribution, we will 
consider two different values for (/34A which we expect to present a reasonable range. 



Next we attempt to relate these unknown initial conditions to the diffractive qq cross section ( 3.13 ) , 
in order to obtain an estimate of their magnitude. The easiest way is a comparison of the two- 
ladder diagrams in AF'^^ with the diffractive cross section (fig. g). We use the coupling ( 2.31| ), 



take the color singlet states of the two ladders, project onto the twist-four term and compute the 
"diffractive" cut o"o between the gluon lines 2 and 3. Equating the result with the cross section 



3.13| we arrive at the equation: 

>4S->4A = ¥'r- (3.21) 



^ 3 jjg 



This result is quite remarkable, since - at first sight - it does not seem to agree with the AGK 
rules. Namely, naively one might have expected that (^943 = ¥54, i.e. the four-gluon contribution 
equals (up to the overall minus sign which we have extracted already in ( p. 17] ) ) the diffractive cross 
section. Let us recapitulate the origin of the l.h.s. of (3.21): from ( 2.16| ) we deduced that QCD 
diagrams for Ty*^* with four gluons in the t-channel have to be decomposed according to their 
symmetry under permutation of color and momenta (the antisymmetric terms have their origin 
in the reggeization of the gluon). The same argument forces us to use the ansatz ( 2.3lD - (|2.33| ), 



i.e. we have the two functions 9943 and (/34a- The AGK rules work for each term separately, but 
the i-channel state with four reggeized gluons couples only to 9943- Finally, when computing the 
two-ladder diagrams in T^*^* , we find for the coupling of two ladders to the proton the combination 
on the l.h.s. of ( 3.21| ) (the fact that the factor in front of 9943 equals 5/4 and not 1 is a result of the 



presence of the second and third term in ( p.32D ). 

There are two important consequences of ( |3.21 ). First, from this single condition we cannot fix 



both initial conditions 9943 and 994a. Secondly, due to the minus sign on the l.h.s., ip^s can be larger 
than the "naive" expectation ip4s = ^4 ■ We propose to proceed as follows. The longitudinal cross 
section for diffractive gg-production allows to estimate 994'^. From our discussion above we know 
that (/?4A is connected with the reggeization of the gluon: in our numerical analysis we will vary 
ip4A in such a way that its correction to the leading-twist structure function ranges between 20 % 
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FIG, 



(a) 



(b) 



(c) 



10: Three sample diagrams which contribute to D4 at order a^ ((a) and (b)) and at order af (c). 
Diagrams (a) and (c) enter also into difFractive qq^-production at leading and next-to-leading 
order, resp. For the difFractive process the two t-channel gluons must certainly be in the color 
singlet, while there are also contributions from other color states which enter into D4. 



and 100 % (see below). Together with the variation of 994A we also find a change in the strength of 
(p4S and AF , the most interesting higher- twist contribution in the small x region. 

3.5 The Sign Structure 

It may be helpful to briefly recapitulate the sequence of our arguments, in particular the sign 
structure. The easiest way is to start with the twist-four corrections to the diffractive cross sections; 
their signs are summarized in tables y and y. Starting with gg-production, in the transverse case 
twist-four is a correction to leading-twist and has a negative sign. The longitudinal cross section 
starts with twist-four and is positive. Next diffractive Q'^g'-production: all signs change, compared 
to gg-production. This is due to the vertex ( 3.15| ) which has an intrinsic minus sign. 



In the next step we have argued that when generalizing to an arbitrary number of gluons in the 
diffractive final state and allowing for evolution in the unintegrated gluon structure function we 
generate all diagrams contained in AF'^^': although there is no direct correspondence between the 
splitting AF^'^' = AF + AF and the different diffractive final states. As a memo we nevertheless 
can say that AF^ "is related to" (i.e. gets the same sign as) the cross section for gg-production, 
and a similar correspondence holds for AF^ and diffractive qqg-, qqgg-, ■ ■ ■ ,-production. 

In the last step we have to remember that it is not only the diffractive states which contribute to 
the structure functions. Making use of the AGK rules we conclude that taking into account all 
other energy discontinuities simply means changing the sign of the (total) diffractive contributions. 
The results are summarized in table Q. 



4 Numerical Results 

4.1 General Procedure and Determination of Initial Parton Distributions 

Using the framework outlined above we have performed a numerical analysis. The first step to 
be discussed is how we fix the free parameters in the various initial distributions (p. From our 
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FIG. 11: The AGK rules: The sum over the three cuts, which is necessary for the computation of the 
inclusive structure function, gives a contribution that is equal to minus one times the term where 
only the diffractive cut is applied. 



discussion it should have become clear that we have to determine four initial distributions, viz. 
1^2 for AF''^^, 973 for the three-gluon correction to /S.F^'^' , and ^p^ which couples to both AF^ and 
AF . Moreover, (^94 consists of a symmetric part and a part with mixed symmetry, cf. eqn. ( 2.3l| ). 
All initial distributions depend on w, and a pole at some positive w-value corresponds to a rising 
(in 1/x) initial condition. 

Let us begin with a detailed discussion of the four-gluon case. As indicated before, we relate (p^ 
directly to the longitudinal diffractive qq cross section (eqs. ( |3.13 ) and ( p.21D ), and the latter can 
be estimated with the help of the measured cross section of longitudinal vector production if we 
sum over all vector mesons in the final state |25]. Let us assume that - after integration over t, 
p1 and M^, and including diffractive dissociation of the proton - this cross section reaches a few 
per cent of the inclusive (transverse) leading- twist cross section aX x = 10~^, Q^ = lOGeV^. Using 
this estimate we fix the initial condition ipf and, with the help of ( 3.21 ), also the combination of 



the initial conditions ip4s and 994A, which enter into the four-gluon amplitude. Since this procedure 
fixes only the difference 5c/94s — 3c/94a (cf. ( p.21| )) we still have freedom for determining the ratio of 
C/94S and (/94A- In a perturbative toy model, namely deep inelastic scattering off a virtual photon, 
(/74A is linked to the leading-twist amplitude: it provides the t-channel four-gluon state inside the 
BFKL amplitude (cf. the discussion in section |3.4D and, hence, can be derived from the size of the 
leading-twist contribution. However, when replacing the (target) virtual photon by the proton, 
this connection becomes somewhat uncertain, and we, therefore, consider different values for the 
ratio ifAh/^AS- Starting with the toy model, we find it reasonable to expect that C/94A should not be 
smaller than 9343. We, therefore, set (^4a = A(/J4s and consider the range 1.0 < A < 1.5. Eqn. ( 3.21 ) 
indicates that this seemingly small variation will nevertheless have a large effect: if A reaches the 
vicinity of 5/3 ~ 1.6, (/54s increases strongly (for fixed y'f )• This enhances the contribution of 
AF^, whereas AF^ changes much less, since it depends on both, (/94s and C/94A (cf. ( p.42| )). Since 
AF and AF enter with different signs, a variation of A between 1 and 1.5 will have a large effect 
on the total twist-four contribution. 

Having fixed the initial conditions (/74 we can, in principle, predict other higher-twist corrections 
to the diffractive cross section. In a future analysis of diffractive jet data at HERA it may be 
possible to measure the twist-four qq and qqg contributions directly; this will provide an important 
consistency check of this strategy. As we have said before, because of the cancellations between the 
qq and qqg parts the total higher-twist contribution to Fg may be smaller than the naive estimate 
based on longitudinal qq production. 

As to the other two initial conditions, 932 and 993, we have not been able to relate them to any direct 
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long, 
trans. 



pT=2 ^p{2) ^^(3) ^^R ^^l 



+ - + - + 

+ + - + - 



TABLE 4: The various contributions to F2 enter with different signs. 



measurement. Starting with c/92, for simpUcity we assume that at the scale Qq the initial distribution 
for twist-four is the same as for the leading-twist case. It follows that the transverse contribution 
is positive, whereas the longitudinal twist-four part of the two-gluon amplitude contributes with a 
negative sign to F2, since 62 in sqn. ( |2.14| ) is negative. For the three-gluon amplitude we make an 
even more ad hoc choice and demand that at x = 10~^, Q^ = 1 GeV^ its absolute magnitude in 
AFf equals the mean value of AF^'^> and AF . Our numerical analysis shows that the choice of 
the point, x = 10~^, Q^ = 1 GeV^, hardly influences our results, since AF^'^\ AF^^^ and AF^ have 
nearly the same shape (up to factors y^aJC^ and as(Q^))- In principle, for both, (/?2^^ and the 
three-gluon term c/33, not only the absolute magnitude but also the overall signs are undetermined. 
We have fixed these signs with the help of the following considerations. First, in AF'^' there is 
no reason why there should be a sign change when going from leading-twist to twist-four. As an 
example, in our toy model described above, deep inelastic scattering of two virtual photons with 
virtualities Q^ <C Q2' both leading- and higher-twist have the same sign. Next, as we have argued 
above, AF^^' and AF^ are supposed to be related to the NLO and NNLO corrections to BFKL. 
Our choice of signs is consistent with the expectation that these corrections have alternating signs, 
i.e. AF^^> is negative with respect to AF^"^' , and AF is positive again. The sign of AF , on 
the other hand, has to be consistent with ( p.42| ) and the AGK cutting rules, i.e. posititive for the 
transverse part. All in all we, therefore, believe that our choice of signs is on a rather safe ground. 
Ta ble ^ summarizes the sign structure of the various leading- and next-to- leading twist parts in F2 , 

cf-dd). 

As one of the results of our analysis we observe that AF^ and AF^ enter with opposite signs and, 
depending on the choice of parameters, tend to weaken each other: although we have some freedom 
for determining the relative contribution of (/94s and tp^A, both, AF and AF , are sizeable and 
enter with opposite signs. But despite these cancellations the total twist-four contribution is not 
small at Q^ = 1 GeV^. It is negative, and the contribution of the transverse photon is substantially 
larger than that of the longitudinal one. 

In our numerical analysis we have used our formulae in order to perform, in the kinematical regime 
of HERA, a "semiquantitative" numerical study of the potential role of twist-four at small x and 
low Q^. As we have emphasized several times the limitations due to the DLA are severe. On the 
leading-twist level it is known that this approximation, which in the anomalous dimension retains 
only the most singular (near u; = 0) term, is insufficient to describe HERA data at small x and low 
Q^: it is, at least, necessary to keep also the constant term within the gluon anomalous dimension. 
We, therefore, expect that a similar level of accuracy will be needed also for a realistic estimation of 
twist-four corrections. In order to minimize the inaccuracy we present the ratios of the twist-four 
contributions to the DLA leading- twist transverse structure function. In this way we hope that our 
results describe the trend correctly on a semiquantitative level. For leading-twist we are fixing the 
free parameter within the initial distribution by demanding that at x = 10~^, Q^ = lOGeV^ our 
DLA expression for F2 is equal to the corresponding MRRS [M value. 
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FIG. 12: Transverse leading-twist structure function in DLA (solid line) in comparison with the longitudinal 
twist-four contribution to F2 of the diffractive gg-production process as a function of x at (a) 
Q^ = 1 GeV^ and (b) Q^ = 10 GcV^. For this scenario, we have chosen constant initial distributions 
at gg = 0.5GeV2. 
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As outlined above, we are using the diffractive vector production cross section to estimate twist- 
four contributions to F2. To be precise, we assume that the longitudinal twist-four qq cross section 
( 3.13| ) reaches three per cent at Q^ = lOGeV^, x = 10~^ relative to the transverse leading-twist 



contribution. As to the choice of the starting scale Qq, our limitations due to the DLA are causing a 
particular problem. In DLA, the four-gluon amplitude /S.F^ vanishes at Q'^ = Qq, and it needs some 
Q^ evolution before it "forgets" about its starting condition. In an analysis which goes beyond the 
DLA, we would start with some nonzero value. In order to escape from this unrealistic behavior 
near Qq, we have chosen a very small value for the input scale, viz. Qq = 0.5 GeV^, and hope that for 
Q^ values greater than, say, 0.8 GeV^, these unwanted effects are suppressed. For the computation 
of the double Mellin transforms we have used the method outlined in ref. |3^, which proves to be 
very accurate and efficient. The numerical results of our analysis also strongly depend upon the 
x-shape of the initial distributions, both for the leading-twist and for the twist-four correction. We, 
therefore, compare three different scenarios, two with constant initial conditions and with different 
choices for the ratio ^aa/'-Pas^ and one with rising initial conditions. 

4.2 Constant Initial Parton Distribution 

Following a presently popular trend in analysing HERA data, we first assume that the leading-twist 
gluon structure function at the input scale Qq is flat at small x. In two subsections we will consider 
the results for two different choices of the ratio A = v'4a/934S) viz. A = 1 and A = 1.5. First, we 



show in figs. 12a and 12b the x-dependence of the transverse leading-twist structure function in 



comparison with the twist-four contribution due to diffractive longitudinal (/^-production. In the 
lower figure the point x = 10~^ is our point of comparison: the higher-twist term has reached 3% 







of F2. The comparison of the two figures shows the expected behavior: higher- twist becomes more 
important at low Q^. If there would be only this twist-four contribution, we would have a positive 
20% correction at Q^ = IGeV^. 

4.2.1 Weak Anti-Symmetric Initial Distribution (A = 1) 

The picture changes quite substantially if we include the other twist-four contributions. We start 
with a discussion of the longitudinal part. Fig. |l^a shows how this twist-four part is composed 
of the four contributions (the twist-four part of the two-gluon amplitude AF'^', AF^'^^ and the 
four-gluon contributions, AF and AF ) at Q^ = IGeV^. In order to simplify comparisons, we 
have normalized all twist-four contributions to the corresponding DLA transverse leading-twist 
values. As can be seen, two of the four contributions (AF^^^ and /S.F ) enter with positive and the 
remaining two with negative signs, so that the sum (solid line) gives a very small contribution. Let 
us take a closer look at the cancellations, in particular those between AF and AF . As we argued 
before, after we have set C/94A = 9743, the relative sign as well as the relative magnitude are fixed: 
in fig.|l3|a (at Q^ = 1 GeV^) AF^ and AF^ almost compensate each other. Due to (|2.45[) AF^ will 



become dominating at very small x. In our calculations this region is not yet reached. 

The sign pattern of the individual contributions does not change when we go to higher Q^-values 
(fig. [13b), although the total sum, now, becomes positive in the full x-interval. A look at the 
absolute values, which are not shown here, tells us that the twist-four contributions have actually 
increased compared to the case at Q^ = 1 GeV^ - the 1/Q^ suppression sets in rather late. But 
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FIG. 13: Longitudinal part of the twist-four contributions (solid line) and the four terms they are composed 
of as a function of x at (a) Q^ — 1 GeV^ and (b) Q^ ~ lOGeV^. Normalization with transverse 
leading-twist values. 
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FIG. 14: Same as in fig. O but, now, for the transverse twist-four contributions. 
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since leading-twist has increased as well, the ratio of the two decreases. Moreover, since AF is the 
twist-four contribution with the strongest evolution (cf. ( |2.45| )) the total twist-four contribution is 
positive and reaches up to 1.5% of leading-twist (fig.|T^). As to the other pair of contributions 
(/\^{2) ^jjjj AF^^'), they also come with opposite signs and, approximately, the same strength. 
But this clearly depends very much on our choice of the initial conditions. For example, a smaller 



choice of (^3, the initial distribution of the three-gluon amplitude, which, as outlined in section 4.1 



cannot be related to any experimental measurement, results in significant smaller values of AFi. 

Next we take a look at the transverse contributions (fig.|l^). As always in our discussion, this 
contribution represents the DLA, and as we discussed before, it is suppressed by one power of 
Iu^Q'^/Qq) in comparison to the longitudinal (twist-four) photon. Nevertheless, our results for 
AFf represent the leading contributions for the transverse structure function. We plot our results 
in the same way as we did for the longitudinal case, and the x distributions for the two values 



Q^ = 1 GeV^ and Q^ = lOGeV^ are shown in figs. 14a and b, resp. First, one notices that for 
each separate twist-four contribution (again normalized to the leading-twist term) the transition 
from Q^ = 1 GeV^ to Q^ = 10 GeV'^ results in a much stronger decrease than in the longitudinal 
case. This is certainly due to the absence of the InQ^-factor, which makes the Q^-evolution weaker. 
Altogether, at Q^ = 1 GeV^ the transverse twist-four terms, which enter into the correction to Ft 
are of much larger size than in the longitudinal case and the sum of all four terms, now, gives a 
negative contribution. 

In our scenario the transverse twist-four contribution (fig.|lja) gives a negative correction which at 
Q^ = 1 GeV^ and x = 10~^ reaches 15% of the leading- twist. The main reason for this is that AF 
is the dominating contribution. If we would go to even smaller values of x the dominance of AF^ 
would become even more pronounced but one has to keep in mind that in this case even higher- 
twist contributions (twist-six, etc.) will, very likely, become important, as well. At Q^ = lOGeV'^ 
(fig. P^b) the transverse twist-four contribution has gone down to less than 4%. Compared to the 
longitudinal case, the twist-four corrections of the transverse photon clearly dominates, and we 
expect that this remains true even if we go beyond the DLA. We can use our numerical results to 



determine the contibution of ip^A within AF^: using ( p. 42 ) and the numerical values in fig. 14 we 



deduce that (for the transverse photon at Q^ = 1 GeV^) it reaches about 20 % of leading- twist. We 
believe that this is a rather low value. In order to further decrease (p4A, we would have to choose 
A < 1. This seems to be rather unlikely if we can believe the results of our toy model. 



Finally, figs. 15a and b show the Q^-dependence of the longitudinal and transverse twist-four contri- 
butions at X = 10"'^. As a general feature, we observe the interplay between the l/Q^ suppression 
factor and the rise in the Q^-evolution. In both figures, the change of the sign in the sum of all 
contributions is due to the increase of AF . Since we have performed our computations in DLA, the 
various twist-four contributions either vanish for Q^ -^ Qq or become a small constant. Moreover, 
AF^ goes much faster to zero than all other contributions. As discussed before, this is a special 
problem of the DLA approximation of AF . If we would take into account also sub-leading terms, 
this decrease of AF^ for low Q^ would be less distinct. It is, therefore, likely that the dominance 
of AF in the twist-four corrections to F2 is, in fact, stronger than found in our DLA analysis. 

In summary, in this scenario the transverse twist-four contribution gives a negative correction of up 
to fifteen per cent of leading-twist at Q^ = 1 GeV^, while the longitudinal twist-four contribution 
remains much smaller. In both cases we find cancellations between the various twist-four contribu- 
tions which enter with different signs. In particular, AF contributes with a rather large amount 
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FIG. 15: Q^-dependence of the (a) longitudinal and (b) transverse twist-four contributions (solid line) and 
the four terms they are composed of at a; = 10^"^. Normalization with transverse leading- twist 
values. 
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compared to AF . In retrospect it seems justified that we have taken into account also the twist-four 
contributions of AF^"^' and AF^^K If we argue that AF^^' should be related to the next-to-leading 
order corrections of the BFKL amplitude, AF^ would correspond even to NNLO corrections, while 
AF^ is the first term which contributes with a stronger power-behavior (cf. eqn. ( |2.45| )). 

4.2.2 Strong Anti-Symmetric Initial Distribution (A = 1.5) 

In this subsection we choose a stronger value for the initial condition (f^A, viz. (p4^x = 1.59343. 
Since iff remains unchanged, our new scenario results in larger values for both 9343 and 934A, 
cf. ( 3.2IJ ). As we have outlined above, AF^ couples only to 9343, while AF^ couples to both 9343 



and 934A. From (2.42) we expect that an increase in 934A leads to a larger value of AF , while 



AF decreases. This is confirmed by the numerical values shown in fig. 16 for the longitudinal and 



transverse contributions at Q^ = 1 GeV^. Both, in the longitudinal and in the transverse case, AF^ 
now dominates. The negative transverse twist-four correction has now gone up to more than 100% 
(at X = 10^'*) of the leading-twist contribution, whereas the positive longitudinal term reaches 
15%. Certainly, if this scenario is realistic, our consideration of twist-four is not sufficient. As soon 
as twist-four reaches, say, fifty per cent of leading- twist, we have to consider also contributions 
with twist larger than four. Using, again, ( p. 42] ) we deduce that in this scenario (for the transverse 
photon at Q^ = IGeV^) the contribution due to 934A reaches about 130% of leading- twist, i.e. 
A = 1.5 is likely to present an upper limit of 934A. 

In both scenarios, considered so far, we find sizeable negative twist-four corrections at Q^ near 
1 GeV'^ : this confirms the expectation that the transition from the parton picture to Regge physics 
at Q^ = is indeed accompanied by sizable higher-twist. In particular, twist-four is negative at 
small X. 

We believe that the largest uncertainties in this analysis are due to the unknown initial condition 
934A. As we have demonstrated, a rather small variation in the ratio 934A/9343 has a large effect on 
the relative contributions of AF^ and AF^, and, therefore, on the total twist-four contributions. 
It is, therefore, important to find further restrictions on the size of (/34A. Other sources of possible 
errors are our choice of the size of AF^^' , and the double logarithmic approximation. Moreover, as 
outlined above, it seems possible that AF^ at very low Q^ would even be larger if one would take 
into account also sub- leading terms. 

4.3 Increasing Initial Parton Distribution 

As an alternative, we have considered a scenario with an initial gluon distribution which increases 
at small x: xg{x,Qfj) ~ (l/x)*^"^. For the twist-four part of the two-gluon amplitude and for 
the three-gluon contribution we assume the same power-behavior of the initial contribution as in 
the leading-twist case, and we have set ujq = 0.2. As to the AF^'^' contributions, the situation 
is slightly more subtle. Since the four-gluon operator in AF^ is closely related to the double- 
exchange of DGLAP ladders, it seems reasonable to assume that the initial distribution of the 
four-gluon amplitude behaves like (l/x)^'^", i.e. the functions 934(0;) in ( p. 31 ) are assumed to have 



a pole at 2ujq. For AF there is an ambiguity: on the one hand, the lowest order diagram is the 
four-gluon exchange, i.e. the evolution should start with 934 which rises as ~ (l/x)^'^°. On the 
other hand, if we consider AF^ as a NNLO correction to the two-gluon operator (|2.2| ) , one might 
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FIG. 16: (a) longitudinal and (b) transverse twist-four contributions as a function of x a.t Q^ = IGeV^. 
For this scenario, we have chosen constant initial distributions with a strong choice of the anti- 
symmetric initial distribution that enters into AF^'^' . The two diagrams, shown here, correspond 
to diagrams |l3| a and |l4| a in the previous scenario. Normalization with transverse leading-twist 
values. 
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FIG. 17: (a) longitudinal and (b) transverse twist-four contributions as a function of a; at Q^ = IGeV^. 
For this scenario, we have chosen increasing initial distributions (wo = 0.2) at Qq — 0.8 GeV^. 
Normalization with transverse leading-twist values. 
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argue that we should have the same initial conditions as for the twist-four part of the two-gluon 
amplitude, i.e. ~ {\/x)^° . In our numerical analysis we have chosen the first option. Finally, since 
it seems unrealistic to assume an increasing initial distribution at Qq = 0.5 GeV^, we, now, set 
Qg = 0.8 GeV^. Moreover we return to our initial convention and set (/34s = ^aa- 



Fig. 17 shows the twist-four contributions at Q^ = 1 GeV^ for the longitudinal and transverse 



photon. Compared to our first scenario, the twist-four corrections are, now, much larger. In order 
to understand this growth of the twist-four corrections we compare the first scenario presented in 



figs. 13a and 14a. Due to the increasing initial distribution, all four contributions are now larger 
than in the case of constant initial distributions. In particular, AF'^' and AF^^> are magnified by 
a similar factor, and the approximate cancellation works in the same way as before. For AF^ and 
AF , on the other hand, the approximate balance is distorted: at Q^ = 1 GeV^ we are close to the 
initial scale, and AF^ is still small. Moreover, because of our choice of the 1/x-power AF^ has 
a stronger rise in 1/x than AF'^^ and /S.F^^> and, therefore, dominates. As result, the twist-four 
corrections at Q^ = 1 GeV'^ are dominated by AF^ which is, now, much less balanced than in the 
previous scenario. By choosing in AF a weaker 2;-dependence at the initial scale and by setting 
our initial scale, as before, to Qq = 0-5 GeV^, it may be possible to partly restore the balance and to 
obtain a somewhat smaller twist-four correction. This scenario, therefore, represents a somewhat 
extreme possibility. 

Nevertheless, the trend seen in this third scenario may be characteristic for increasing initial dis- 
tributions. Initial parton distributions which increase at small x tend to increase the twist-four 
corrections at low Q^, in particular through AF^. This pushes the balance between the positive 
AF^ and the negative AF^ into the positive direction, and, as in our case, may even lead to a net 
positive twist-four correction to the leading-twist structure function. If these tentative conclusions 
are correct, we have to interpret the scenario with rising initial conditions as being unrealistic. As 
discussed in the beginning of this paper, twist-four corrections could, in principle, provide an ex- 
planation for the observed change in the small x behavior at low Q^. In particular, with decreasing 
Q^, a negative transverse higher-twist could, because of its stronger rise at small x, overcome the 
l/Q^ suppression and interfere with the rising leading-twist contribution. In this last scenario, 
however, we found that the opposite sign for transverse higher-twist is more likely, and leading- 
and higher- twist are adding up rather than compensating each other at low Q^ . We believe that 
such a scenario is not very likely to be realistic. 



5 Conclusions 

In this paper we have made a first attempt to estimate twist-four corrections to the gluon structure 
function at small x and low Q^. We have collected all presently available analytic information on 
the twist-four gluon evolution equations and defined the terms which need to be included. Up to 
now DLA is the only approximation for which explicit expressions exist. One of the main features 
of gluonic twist-four corrections is the sign structure: we find, both in the transverse and in the 
longitudinal structure functions, terms with opposite signs which tend to compensate each other. 

In the second part we have investigated how DIS diffractive dissociation (in particular the produc- 
tion of longitudinal vector mesons) can be used for constraining the free parameters, i.e. the initial 
conditions. As a main result, the AGK rules provide a bridge between DIS diffractive dissocia- 
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tion and the twist-four part of the four-gluon ampHtude, but there remams a piece in the initial 
condition which cannot be determined in this way. 

In our numerical part we have considered three different scenarios. For two of them we have chosen 
constant (in x) initial parton distributions at the input scale Qq, and we have varied the relative 
strength of pieces in the input distribution. In addition, we have studied also one scenario in 
which the input distribution rises in 1/x. In general we find, both in the transverse and in the 
longitudinal structure functions, significant cancellations between terms with opposite signs which 
tend to weaken or even compensate each other. Moreover, the transverse twist-four corrections are 
larger than the longitudinal ones. For the scenarios with constant initial conditions the sum of 
the transverse twist-four corrections is negative, and at small x grows stronger than leading- twist. 
For Q^ = 1 GeV^, x = 10~^ twist-four ranges between 15% and 130% of leading- twist. Compared 
to the results of [^, our corrections seem to become significant in the right region. For the case 
of increasing initial conditions the different pieces of the transverse twist-four contributions add 
up to a positive higher-twist correction; we interpret this scenario as rather unrealistic. This may 
however change if one could go beyond the DLA. All our numerical results are based upon DLA 
accuracy and, therefore, should be taken as only "semiquantitative" . 

Before firm conclusions can be drawn, several improvements have to be made. First of all, we 
need the complete LO calculations of the evolution kernels of the gluonic operators. Work along 
these lines is in progress [53], but it will take some time before one can start to perform numerical 



calculations. Secondly, as can be seen already from eqn. ( 2.38| ), the amount of computer work 
necessary to calculate the Q^ evolution of the gluon operators is substantially larger than in the 
leading-twist case. Presently, it is not clear whether it is possible to find a closed expression also 
beyond DLA or whether we have to solve evolution equations in x-space. Attempts to solve the 



set of coupled evolution equations (2.1§|) will face difficulties with storage capacities. Finally, as 



discussed above, experimental information on specific final states is needed in order to restrict the 
arbitrariness in the choice of the initial conditions. In this first attempt we have used only the 
diffractive longitudinal vector meson production in order to obtain a first hint on the magnitude 
of twist-four. Besides this, other information of final states might be useful; a promising candidate 



might be a measurement of the two-particle correlation function as a function of rapidity |36|. 
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Appendix 

We, first, give expressions for the vectors in color space \i) in the couphng scheme x = (1234): 
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and eigenstates of the BFKL color factor t2^2 : 
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and Ei = {— 3, — |, — |, 0, 1} for i = {1,8a, 8s, 10 + 10,27}. They are related to vectors in other 
coupling schemes via transitions with the matrices 
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and P = diag(l,-1, 1,-1,1) (A.9) 



in the following way: 
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Next, we give explicit expressions for the components of vector 

x{ujD) = T.{uji>)V (A.ll) 

in eqns. ( |2.3[1| ) and ( p.38| ). The results are the same for the three coupling schemes. Therefore, we 
list the components of x ^^r the three contributing color states and write them as functionals of 
the auxiliary quantities fi(ujD), i = 1 . . .6, and gj{uji'), j = 1 . . .7. They are defined as follows: 

_J_ A_^ _1 f-1 L_J__ZL -1. 

■'^ ~ ~32^^ ^ 40^^ ~ 160^^^ ~ 4 71 ~ 128^^ ~ 16^® ~ 128^^"^ ~ 16 

■^^ 128^^ 80^® 640^^^ 16 78 16^^ 200^® 400^^"^ 20 ^ ' ' 

3 9 21 3 ^ 1 27 27 49 7 

/3 = -TTTTTI + 7:7^78 - 7^7:7:727 - 7:7: J6 = 77:7771 " -77:7778 - 777:777:727 - 7777 



and 



32' 200' 800' 20 727 128' 400' 3200' 80 



91 = ^/i • /s - 3\/3/2 ■ h 94 = h- h- -^fi 

92 = ^/2 • /3 - ^/i • /6 55 = /4 • /6 - 27/| (A.13) 

3\/6 ..\/6,. i^i^l^2 

53 — ~^/l ■ /2 ^/S • J4 56 — 74 • /5 - -^/l 

57 = 9/1 ■f2-h + U-h-k- \fl ■ k - 27 f^ ■ h - \fl ■ h , (A.14) 



where the functions 7j(u;i>) are the diagonal elements of matrix ^{ujD) in eqn. ( 2.28| ), and we have 
used 78 instead of 78g . With these functions the elements of the vector x are given by the following 
expression: 



57 • Xj = Vj97lj + Y.9i ('^S?7i7i + 4i^ljl8 + 4^ilj727 + 4? 71 + 4? 78 + 1^6-127) , (A.15) 



yi \r ■ ' 

i=l 

where j G {1, 83, 27}, vi = 2, vg = 4-v/2, i;27 = 6\/3 (cf. ( |2.2g| )) and the coefficients kJ^ are listed in 
table ^ The pole at u>u = 4(1 + 6)Ncas/TT (eqn. ( ^.39|) ) is due to a zero of function 57, eqn. ( |A.14 ). 

Finally, the color tensor d"''"^'^ in ( |2.41[ ) has the form p: 

_/i'*^A'^'^ _l_ _^cibljlcd _ j:abl Acd 

~ 6 4 4-' -^ 

= ljad^6c ^ \^adl^lhc ^ \ pdl jlbc ^ ^p^-^Q-^ 

6 4 4 

where the Gell-Mann matrices t"" are normalized to 

tr{n''} = -(5'^^ (A.17) 
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i 


1 


2 


3 


4 


5 


6 




3v^ 
32 


16 


3v^ 
16 


1 

64 


1 

8 


27 
64 




9v^ 
20 


7V2 
40 


3^6 
40 


1 
8 


3 

10 


27 
40 


^3J 


213v^ 
160 


81v^ 
80 


51v^ 
80 


27 
64 


27 
40 


189 
320 




3v^ 

4 


2 




1 
4 








6v^ 
5 


6^2 
5 




2 








21v^ 
20 


27v^ 
10 




27 

4 








9^6 
80 


7 
80 


3v^ 
80 


V2 
32 


3^2 
40 


27v^ 
160 


'*2J 


3^6 
10 


3 

~5 


9^3 
25 


4 


9^2 
50 


27^2 
100 


^3J 


51v^ 
80 


27 
80 


99^3 
400 


27v^ 
32 


81v^ 
200 


189^2 
800 






1 
4 


3v^ 
4 




V2 
2 








2 


6v^ 
5 




6^2 
5 




^6i 




27 
4 


21v^ 
20 




27v^ 
10 




^(27) 


71 
160 


9v^ 
80 


17v^ 
80 


3v^ 
64 


3v^ 
40 


21v^ 
320 


^(27) 
'*2J 


17 
20 


3v^ 
40 


33 v^ 
200 


3v^ 

8 


9v^ 
50 


21v^ 
200 


^(27) 
'*3J 


189 
160 


81v^ 
80 


189v^ 
400 


81^3 
64 


81^3 
200 


147v^ 
1600 


^(27) 
'*4i 


1 

4 




V2 
2 






3^3 
4 


^(27) 


2 




6^2 
5 






6^3 
5 


^(27) 


27 
4 




27^2 
10 






21v^ 
20 



TABLE 5: Expansion coefficients for the three different color states in eqn. (A. 15) 
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